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We present the first study of the energy dependence of p t angular correlations inferred from 
event-wise mean transverse momentum (pt) fluctuations in heavy ion collisions. We compare 
our large- acceptance measurements at CM energies ^/sjvjv = 19.6, 62.4, 130 and 200 GeV to 
SPS measurements at 12.3 and 17.3 GeV. p t angular correlation structure suggests that the 
principal source of p t correlations and fluctuations is minijets (minimum-bias parton fragments). 
We observe a dramatic increase in correlations and fluctuations from SPS to RHIC energies, in- 
creasing linearly with In ^/sjvjv from the onset of observable jet-related (p t ) fluctuations near 10 GeV. 
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I. INTRODUCTION 



II. ANALYSIS METHOD 



Theoretical descriptions of heavy ion (HI) collisions 
at RHIC energies predict copious parton (mainly gluon) 
production in the early stages of collisions and subse- 
quent parton rescattering as the principal route to a 
color-deconfined bulk medium, with possible equilibra- 
tion to a quark-gluon plasma [J, [2, EJ El, 111- Parti- 
cle yields, spectra and high-p t correlations from Au-Au 
collisions at 130 and 200 GeV provide tantalizing evi- 
dence that a QCD colored medium is indeed produced at 
RHIC SSBEEEi- What are the properties of that 
medium, and has an equilibrated QGP formed prior to 
hadronic decoupling? A partial answer may emerge by 
searching for evidence of initial-state semi-hard scattered 
partons in the correlations and fluctuations of final-state 
hadrons. In particular, p t fluctuations and correlations 
may provide such evidence. 

Fluctuations (variations about a mean) of event- 
wise mean transverse momentum (p t ) [l2l Il3| within 
momentum-space angular bins of varying size, and corre- 
sponding two-particle pt correlations (variation of a two- 
particle distribution relative to a reference), could pro- 
vide access to early parton scattering and subsequent in- 
medium dissipation inaccessible by other means [Til ]. p t 
correlations may represent such partons in the form of 
local velocity and/or temperature correlations [HI, [lj|. 
Measurements of (pt) fluctuations in Au-Au collisions at 
fixed scale (bin size) at 130 GeV [l2[ and measurements 
of p t angular correlations inferred from (p t ) fluctuation 
scale dependence at 200 GeV pj} indicate that p t corre- 
lations at RHIC are much larger than those at the SPS. 

In this paper we report the first study of the energy 
dependence of pt angular correlations (e.g., structures in 
the event-wise pt distribution on (77, <p) which occur at 
different positions in each HI collision) inferred from ex- 
cess (pt) fluctuations (fluctuations beyond those expected 
for independent particle p t production). We present the 
scale dependence of (p t ) fluctuations within the STAR 
detector acceptance for four RHIC energies and pro- 
vide a basis for interpreting those fluctuations by in- 
verting the fluctuation scale dependence for two energies 
and two centralities to form p t autocorrelations on angle 
space (77, (/>). We obtain the centrality dependence of full- 
acceptance fluctuations at four RHIC energies compared 
to results at two SPS energies, and we determine the 
y/sNN dependence of (pt) fluctuations for full-acceptance 
STAR data as a basis for comparison with extrapolated 
CERES measurements [181 and the pQCD event simula- 
tion Monte Carlo Hijing This analysis is based on 
Au-Au collisions observed with the STAR detector at the 
Relativistic Heavy Ion Collider (RHIC). 



Excess charge-independent (all charged parti- 
cles combined) (pt) fluctuations are measured by 
the difference between the variance of quantity 
{pt(5x) — n(5x)pt}/ \Jn(8x) and the variance tr| t 
of an uncorrelated reference [l2j]. pt(Sx) is the scalar 
sum oi p t in a bin of size 5x (e.g., 5rj or 5<f), n(8x) is the 
number of particles in the bin, and pt is the mean and 
Cp t is the variance of the single-particle pt spectrum for 
all accepted charged particles from all events (pt and cr| t 
then represent independent particle p t production) . The 
variance difference 



&(T pt:n (5x) = {p t (Sx 

= 2a Pt Aa Pt:n (Sx), 



n(6x)pt} 2 /n(6x) — a\ 



(1) 



calculated over a range of bin sizes is the scale-dependent 
(p t ) fluctuation measure for this analysis [l7|. Overlines 
indicate averages over all bins with size Sx in all events. 
The difference factor Aa pt . n is related to the Opt) fluctu- 
ation measure $ pt introduced previously [2(| l2ll . |22| by 
Aa Pt:n = <P Pt (l + <i> Pt /2a Pt ). 

Event-wise fluctuations in the bin contents of a binned 
distribution reflect changes in the structure of that distri- 
bution. Smaller bins are sensitive to more local aspects 
of that structure, and conversely. The scale dependence 
of fluctuations is therefore equivalent in some sense to a 
running integral of the corresponding two-particle distri- 
bution, the integration limit determined by the bin size or 
scale [23I l24j . The integrand is an autocorrelation which 
compares a distribution f(x) to itself. An autocorrela- 
tion is effectively a projection by averaging of product 
distribution f(x\) ■ f(x2) on (x\,X2) onto the difference 
variable a; a = x\ — a; 2 [2q . In this analysis we wish to de- 
termine the average angular correlation structure of the 
event- wise pt distribution on (77, (f>): what are the aspects 
of that distribution which vary event-wise but which nev- 
ertheless have persistence and universality. We do so by 
inverting the scale dependence of excess (pt) fluctuations 
(the integral) to obtain the autocorrelation of the pt dis- 
tribution on (77, 0) [13, H3, HI] • 

Eq. ^ below is an integral equation in discrete form 
which relates variance difference Aa pt . n (8r], S4>) on pseu- 
dorapidity 77 and azimuth angle <p to an autocorrelation 
distribution on (n A ,4>A) (e.g., 77 A = m ~ V2) [13, HI- 
The autocorrelation (cf. Fig. [5] for examples) compactly 
represents two-particle correlations on (77, 4>) in HI colli- 
sions !26j . The 2D discrete integral equation is 



(2) 



. K &p(Pt : n;ke v ,le<i,) 

k,l=l 



Typ r ef(n;ke v ,le<j>) ' 



with kernel K mn .ki = (m — k + l/2)/m ■ (n — I + l/2)/n 
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and fixed microbin sizes e n and for the discrete 
integral. That equation can be inverted (solved for 
the integrand) to obtain autocorrelation density ratio 
Ap{p t : n)/ y/p re f(n) [units (GeV/c) 2 ] as a per-particle 
p t correlation measure on (ryAj^A) fro m (p t) fluctuation 
scale dependence of Ao-p t . n (5r],5(j)) [TEHl]- Ap(p t '■ n ) is 
proportional to the average of (pt—npt) covariances for all 
pairs of bins (^Aj^a) apart, y/ p re f (n) is the geometric 
mean of particle densities in those bins [13] • Density ra- 
tio Ap{p t : n) I \J p re f(n) is thus proportional to normal- 
ized covariance (pt — np t ) a {pt — npt)^/ ' \fn^n^ (averaged 
over certain bin combinations (a, b) to form an autocor- 
relation). The densit y ra tio has the form of Pearson's 
correlation coefficient (2 7} , but with number variances in 
the denominator replaced by Poisson values n a , fib. The 
density ratio is derived and discussed in [H, [2f| [2i| [3(| ■ 



III. DATA 

Data for this analysis were obtained with the STAR de- 
tector [3l[ using a 0.5 T uniform magnetic field parallel 
to the beam axis. Event triggering and charged-particle 
measurements with the Time Projection Chamber (TPC) 
are described in (3l| . Track definitions, tracking efficien- 
cies, quality cuts and primary-particle definition are de- 
scribed in [12|, HH ■ Tracks were accepted with pseudora- 
pidity in the range |7y| < 1, transverse momentum in the 
range pt €E [0.15,2] GeV/c and 2ir azimuth, defining the 
detector acceptance for this analysis. Particle identifica- 
tion was not implemented. Centrality classes (percent- 
ages of the total hadronic cross section) were defined in 
terms of the uncorrected number N of charged particles 
in acceptance \r]\ < 1 according to procedures described 

in dilizl. 

Centrality specified in terms of pathlength v (esti- 
mating the mean number of nucleons encountered by 
a participant nucleon) is based on the relationship of 
N to minimum- bias distribution endpoints N p and Nq. 
N p , the lower half-maximum point of the minimum- 
bias distribution plotted as do/dN 1 / 4 , estimates the 
uncorrected mean multiplicity for non-single-diffractive 
nucleon- nucleon collisions in the same acceptance. No, 
the upper half-maximum point of do / dN 1 / 4 , estimates 
the value of N corresponding to the maximum number 
of participant nucleons N par t, max and impact parame- 
ter b = 0. The relation between the fractional cross 
section o/oq and multiplicity N is approximated by ex- 
pression 1 - o/o = (N 1 / 4 - iV p 1/4 )/(7V V4 - 7V p 1/4 ) [H, 
accurate to ~ 3% over the entire centrality range (exlud- 
ing fluctuations near the endpoints). Mean participant 
path length v w 2Nun/N part [3J| is then given by 
v = {1+2.23 (1-ct/ct )} 6 /{(1+2.72 (1-o/o )} 4 to about 
2%, based on a Monte Carlo Glauber simulation [33| . Co- 
efficients 2.23 and 2.72 apply to 200 GeV collisions and 
vary slowly with ln(y / s]vjv), resulting in a few-percent 
shift in v max with collision energy over the energy range 




FIG. 1: Per-particle (p t ) fluctuation scale dependence for 
19.6, 62.4, 130 and 200 GeV Au-Au collisions, and for the 
top 0-20, 0-5, 0-15 and 0-5% of the total hadronic cross sec- 
tion respectively. 



62 - 200 GeV. Those expressions determined the values 
of v used in Fig. [3] (right panel). 



IV. FLUCTUATIONS AND CORRELATIONS 

Fig. Q] compares the scale dependence of variance dif- 
ference Ao-p t . n (Sr],5(f)) in Eq. Q] for central Au-Au col- 
lisions and four collision energies: ^/snn = 19.6, 62.4, 
130 and 200 GeV. The increase of fluctuation amplitudes 
between 19.6 and 200 GeV in Fig. [1] is a factor four, 
establishing that the (pt) variance difference is strongly 
energy dependent. However, (p t ) fluctuations are diffi- 
cult to interpret, whereas the corresponding p t angular 
autocorrelations obtained by inverting fluctuation scale 
dependence clearly indicate the underlying dynamics. 

Fig- HI shows pt angular autocorrelations (by construc- 
tion symmetric about ?7a, 0a = 0) of density ratio 
Ap/ y 'p re f inferred from (p t ) fluctuation scale depen- 
dence as in Fig. [1] by inverting Eq. ([2]) [23| . The plots 
in Figs. [1] and [2] contain equivalent information in differ- 
ent forms, consistent with Eq. ^j. One can observe the 
equivalents of the elliptic flow sinusoids in Fig. [2] (right 
panels) along the upper-right edges of the plots in Fig. Q] 
Note that the elliptic flow sinusoid amplitudes in this 
autocorrelation representation using per-particle correla- 
tion measure Ap/ y/p re f are negligible for peripheral col- 
lisions and increase with increasing centrality toward a 
maximum for mid-central collisions. That behavior con- 
trasts with the nearly opposite trend observed with con- 
ventional per-pair measure i>2- The top two panels of 
Fig. [U represent 130 GeV and the bottom two panels 62.4 
GeV Au-Au collisions. The left panels represent compa- 
rable peripheral collisions and the right panels compara- 
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FIG. 2: pt autocorrelations for 130 GeV (upper) and 62.4 
GeV (lower), and for comparable peripheral (left) and cen- 
tral (right) Au-Au collisions. These autocorrelations contain 
same-side (|0a| < tt/2) and away-side (\</>a\ > 7r/2) struc- 
tures. 



circumstances (e.g., detector acceptances) may result in 
apparently conflicting fluctuation measurements. How- 
ever, per-particle fluctuation measurements such as those 
presented here are exactly comparable at the same bin 
size or scale [integration limits of Eq. @], indepen- 
dent of detector geometry and other experimental de- 
tails, because they integrate the underlying autocorrela- 
tions which are detector-independent distributions. (p t ) 
fluctuations have been measured by several collabora- 
tions [H [H, EE EEJp, HI • PHENIX measurements at 
130 and 200 GeV~[l3l] are compatible with STAR mea- 
surements at equivalent acceptances (scales). We wish 
to determine the energy dependence of pt angular cor- 
relations over the largest energy interval possible, from 
200 GeV at RHIC down to the lowest SPS energies. How- 
ever, the CERES SPS measurements are restricted to the 
77 scale dependence of (pt) fluctuations at full azimuth. 
We make the most differential comparisons possible with 
that limited information, given that fluctuation scale de- 
pendence is the running integral of the underlying au- 
tocorrelation. First, we make a detailed comparison of 
STAR and CERES fluctuation scale dependence at full 
azimuth and equivalent pseudorapidity scales. Then we 
examine changes in centrality dependence with energy. 



ble central collisions. Autocorrelations for 200 GeV are 
presented in and the available 19.6 GeV data do not 
have sufficient statistics for a satisfactory fluctuation in- 
version. The autocorrelations in Fig. [2] have same-side 
(|0a| < 7r /2) and away-side (|0a| > 7r /2) components. 

We interpret the same-side peak and away-side ridge 
which dominate peripheral collisions as consistent with 
minijets as the source mechanism. The shapes are simi- 
lar to jet correlations observed in p-p collisions [30| . and 
in Au-Au collisions of all centralities as modeled by Hi- 
jing [23 | . Those structures are strongly modified with 
increasing A-A centrality, but in a continuous manner 
which suggests that the minijet interpretation is also ap- 
propriate in central heavy ion collisions. By 'minijet' we 
refer to correlated hadrons from initial-state semi-hard 
parton scattering and subsequent fragmentation in which 
no leading or trigger particle is required, i.e., fragments 
from minimum-bias partons with no analysis restriction 
placed on the parton momentum spectrum. We expect 
minimum-bias partons to be dominated by the low-Q 2 
'minijets' described by theory [2, H, H, The ampli- 
tudes of the peripheral same-side peaks are notably sim- 
ilar for 62.4 and 130 GeV, whereas for central collisions 
the same-side peak amplitude increases strongly with en- 
ergy, and the peak is significantly broadened in the r/A 
direction. 



A. Scale dependence 

Fig. [3] (left panel) shows CERES $ Pt values for Pb-Au 
collisions at several pseudorapidity scales and 27r azimuth 
acceptance for ^/snn — 12.3 and 17.3 GeV IS]. Also 
plotted are comparable STAR measurements of Ao~ Pt:n . 
The CERES data rise rapidly to about 3 MeV/c within 
Sr/ < 0.2 (see inset), with a slower linear rise thereafter. 
STAR 19.6 GeV data show similar behavior, albeit with 
somewhat larger magnitudes over a larger r\ acceptance. 
The higher-energy STAR data are qualitatively larger in 
magnitude. We attribute the rapid rise of CERES data 
(~ 0.003 GeV/c) in Sn < 0.2 to quantum (HBT) and 
Coulomb correlations (resonance decays make a negligi- 
ble contribution to (p t ) fluctuations) and designate those 
contributions as small-scale correlations (SSC). The com- 
plementary region of Sn then represents large-scale cor- 
relations (LSC). We conclude (see below) that the LSC 
component is dominated by parton fragments (including 
the away-side ?7A-independent azimuth peak) , and possi- 
bly global temperature fluctuations [with corresponding 
autocorrelation uniform on (t^a^a)]- 

B. Centrality dependence 



V. EXPERIMENT COMPARISONS 

The autocorrelations in Fig. [2] represent all the angu- 
lar correlation information obtainable from correspond- 
ing fluctuation measurements. Different experimental 



Fig. [3] (right panel) shows the centrality dependence of 
charge-independent (pt) fluctuations (elliptic flow does 
not contribute to (p t ) fluctuations integrated over a 2tt 
azimuth acceptance). Ao~ Pt:n was measured at the STAR 
acceptance scale and at four collision energies and cor- 
rected for tracking inefficiency and background contam- 
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FIG. 3: Left panel: Per-particle fluctuation dependence on 
pseudorapidity scale Srj (in 2n azimuth) in central collisions 
(c/. Fig. [T] caption for STAR centralities) . STAR measure- 
ments are solid symbols, CERES measurements [18J are open 
symbols. The inset shows details at small Srj. Right panel: 
Centrality dependence of (pt) fluctuations in the STAR ac- 
ceptance for four energies, v is the mean participant path 
length 34]. The vertical line at right estimates v for b = 0. 
The upper hatched band estimates the uncertainty in v for 
130 GeV data. There is an overall 14% systematic error in 
the corrected amplitudes. SPS measurements of "J> pt at 12.3 
and 17.3 GeV (the lower-right hatched region, with errors and 
centrality range) are included for comparison. Curves guide 
the eye. 



VI. ERRORS 

Statistical errors for (p t ) fluctuation data in Fig. Q] are 
0.005, 0.001, 0.003 and 0.0015 (GeV/c) 2 respectively for 
19.6, 62.4, 130 and 200 GeV (statistical variations at dif- 
ferent scale points are correlated by the integral nature 
of fluctuation scale dependence). Errors for fluctuation 
measurements in Fig. [3] (right panel) are indicated by er- 
ror bars. The upper hatched band in Fig. [3] (right panel) 
estimates the uncertainty in v for 130 GeV data and pro- 
vides an upper limit for 62.4 and 200 GeV v uncertainties 
(more typically < 3%). The uncertainty in v for the 20 
GeV data is approximately 0.3, as shown by horizontal 
errors in Fig. [3] (right panel) . Systematic corrections to 
fluctuation amplitudes for tracking inefficiency and back- 
grounds vary over 15-22% and 21-35% ranges for 130 and 
200 GeV respectively. The overall systematic uncertainty 
for corrected fluctuation amplitudes is 14%. Autocorre- 
lation errors have two components: statistical fluctua- 
tions which survive smoothing and systematic error due 
to smoothing distortion. Statistical errors for the auto- 
correlations, estimated by inverting the error estimate 
for Aa^ t:n , are less than 0.0001 and 0.0003 (GeV/c) 2 for 
autocorrelations at 62.4 and 130 GeV respectively (bin 
errors are correlated). Smoothing distortions, estimated 
by passing data through inversion twice, are less than 5% 
of the range of autocorrelation values in each panel. 



ination as in jl2j. The vertical scale in this figure rep- 
resents r.m.s. fluctuations measured by Aer pt:n (STAR). 
Corresponding $ Pt (CERES) values are numerically the 
same within 1.5% (less than the relative errors) within 
the range of the CERES data [HI , and are therefore plot- 
ted without correction on the same scale. 

The fluctuation amplitudes in Fig. [3] (right panel) vary 
strongly with collision centrality (consistent with pjj]) 
and energy. The observed trends are fully consistent 
with p t angular autocorrelations reported in (specif- 
ically the near-side minijet peak amplitude). Although 
there is a trend of monotonic increase with energy for 
the more central collisions, there is an interesting satura- 
tion of y/SNN dependence for peripheral collisions (v < 
2.5), consistent with the similarity between left panels in 
Fig. [21 Also included in this panel is a summary (lower- 
right hatched box) of $ Pt measurements for Pb-Au col- 
lisions at 12.3 and 17.3 GeV extrapolated to the STAR 
rj acceptance (see left panel) for comparison. The lower- 
left hatched box at v = 1 represents a $ Pt measurement 
of 2.2 ± 1.5 (sys) MeV/c for p-p collisions at 17.3 GeV 
(in the forward rapidity acceptance € [1.1,2.6]) [22| . 
We attribute no special significance to the apparent gap 
between 19.6 and 62.4 GeV data in Fig. 3. In Fig. g] 
the data are consistent with the simple logarithmic trend 
ln(y'sjvAr). The gaps in Fig. [3] arise from a conspiracy of 
currently available energies and the logarithmic trend. 



VII. ENERGY DEPENDENCE 

As noted previously, although the energy dependence 
of p t angular autocorrelations has been established for 
STAR data by inversion of (p t ) fluctuations (e.g., Fig. [2] 
and [13]) Pt autocorrelations have not been measured at 
SPS energies. Therefore, we infer the energy dependence 
indirectly using (pt) fluctuation measurements as proxies. 
The energy dependence of (pt) fluctuations for STAR and 
CERES data is summarized in Fig. Fluctuations mea- 
sured by the most-central, full- acceptance STAR Aa Pt:n , 
and CERES $ Pt values linearly extrapolated to 6r) = 2 
(CERES data are linear on Srj in [0.3,0.8], cf. Fig. [3]- left 
panel) , are plotted in the left panel vs ^/snn, with (solid 
points) and without (open points) SSC correction. As lin- 
ear extrapolations the CERES points are actually upper 
limits, since an NA49 measurement at 17.3 GeV in the 
pion rapidity interval [1.1,2.6] in the CM gave an upper 
limit for $ Pt of 1.6 MeV/c for central collisions [Hf. A 
subsequent measurement in the same rapidity interval re- 
vealed nonzero results for more peripheral collisions [22| . 
The SSC correction consists of subtracting 0.003 GeV/c 
from CERES $ Pt and STAR Acr pt: „ values (the contri- 
bution to the integral of Eq. for Srj < 0.2). Given 
the uncertainly in the SSC correction and the CERES 
extrapolations the open points at and below 19.6 GeV 
represent upper limits on a minijet contribution. The 
solid curve, proportional to ln{^/sAr at/10}, summarizes 
the trend of the data above 10 GeV. We observe a dra- 
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FIG. 4: Left panel: y^ijvjv dependence of (p t ) fluctuations for 
central collisions and STAR full acceptance. CERES fluctu- 
ation data (12.3 and 17.3 GeV) were linearly extrapolated to 
the STAR rj acceptance (cf Fig.[5J. The curve is proportional 
to ln{y / ijvjv/10}. Right panel: Hijing-1.37 (pt) fluctuations 
for central collisions and the STAR acceptance at three ener- 
gies, for quench-on, quench-off and jets-off collisions. 



matic increase of (pt) fluctuations with beam energy for 
corrected and uncorrected data from an apparent onset 
of minijet-related correlation structure near 10 GeV. 

The fluctuation measure £ 



Pt - y Acr p t :n/ n P? (quanti- 
ties in the radicand are as defined in this paper) was em- 
ployed in [lU to suggest that the energy dependence of 
(p t ) fluctuations is negligible, dominated by global tem- 
perature fluctuations with o~t/Tq ~ 1%, unchanging from 
SPS to RHIC [H, [Hj]. £ Pt is based on four assumptions: 
1) each collision event is thermalized with temperature 
T, 2) an event ensemble has nonzero temperature vari- 
ance er|, about ensemble mean To, 3) Acr pt . n /n estimates 
dip and 4) p t estimates To. £ Pt should therefore estimate 
ax /To. All four of those assumptions are falsified by the 
results reported in this paper. The event-wise thermal- 
ization scenario implied by assumptions 1) and 2) is falsi- 
fied by the strong minijet pt correlations shown in Fig. [2] 
and in more detail in [1 71 ] . The observed pt correlations 
are almost three times larger than those predicted for no 
thermalization of jets (quench-off Hijing) in 24]. The 
separate elements in £ Pt do not estimate temperature- 
related quantities as implied by assumptions 3) and 4). 
We learn from Fig. [5] and more extensive results in [17| 
that A(jp t:n is dominated by jet correlations when eval- 
uated at 8(f) = 2ir (the elliptic flow contribution then 
integrates to zero). The issue of parton scattering and 
incomplete equilibration is discussed further in Sec lVIIll 
Each of h, Acr pt . n and p t varies strongly with collision 
centrality and energy. Measured angular p t autocorre- 
lations combined with two-component trends for n and 
Pt strongly suggest that those variations are dominated 
by incompletely-equilibrated semi-hard parton scatter- 
ing. If the observed hard-scattering contribution were 
more strongly equilibrated S P4 might even fall sharply 
with increasing energy. The algebraic combination £ Pt as 
measured happens to nearly cancel the ^/s n Ar-dependent 
hard-scattering trends of the individual factors, but the 
significance is unclear. A clearer picture emerges when 



the quantities are studied separately, as in the present 
study. S P( is dominated by reference factor 1 / -J N part 
and mixes soft (HBT, Coulomb) and hard (parton frag- 
ment) contributions (SSC and LSC components respec- 
tively) through a running average on scale Sr/. £ Pt is thus 
by construction insensitive to the observed energy depen- 
dence of parton scattering and fragmentation which dom- 
inates (pt) fluctuations at RHIC and which motivated 
this study. 

An analysis using a variant of £ Pt denoted 
\J {o~Pt,i ■ Spt,j) I ((pt)) was reported in (36|. However, 
whereas CERES' E Pf is constructed with the same under- 
lying statistical quantity Aer pt . n used in the present anal- 
ysis, and is therefore directly comparable, the quantity 
\J (Sp t ,i ■ 5p t j)/((pt)) is not. In particular, (5p t .i ■ 5p t .j) 
includes in its denominator the random variable n(n — 1) 
(particle multiplicity n varies randomly from event to 
event within some limits) [12]. For that reason it can 
produce results significantly inconsistent with $ Pt and 
AtTp t . ra when bin multiplicities are small. 

For example, in Fig. [2] (lower-left panel) we show re- 
sults for 80-90% central Au-Au collisions. The corre- 
sponding structure for the Hijing Monte Carlo is very 
similar: a same-side minijet peak and an away-side ridge 
from back- to-back jets [24]]. However, the minijet peak 
from n (5pt,i • $Pt,j) applied to the same Hijing data is 
about 3x smaller, and there is additional (mainly nega- 
tive) structure of comparable magnitude not present in 
the A(7p t . n result. Deviations in Hijing data become sig- 
nificant for centralities more peripheral than 50%. Thus, 
(Spt.i-Sptj) cannot be used to study the scale dependence 
of pt fluctuations, where the average multiplicity in a bin 
can be as small as 1. 



Analysis of (pt) fluctuations for Hijing-1.37 central col- 
lisions in the STAR acceptance at 200 GeV [H [II 
was extended down to 62.4 and 19.6 GeV for this 
energy-dependence study. Values of A<r Pt: „ for quench- 
on, quench-off and jets-off central Hijing collisions (jet 
quenching modeled by pQCD gluon bremsstrahlung) are 
shown in Fig. 2] (right panel). The variation with energy 
of Hijing (p t ) fluctuations is small for jets-on (quench-off 
and quench-on) collisions. However, the corresponding 
Hijing pt autocorrelations (cf. [24[ for examples at 200 
GeV) reveal that the quench-off, near-side minijet peak 
amplitude falls by about lOx from 200 to 19.6 GeV, con- 
sistent with the variation of (pt) fluctuations from data 
shown in the left panel. The source of the inconsistency 
is apparently the Hijing string fragmentation model, pt 
correlations in Hijing below pt — 0.5 GcV/c from string 
fragmentation persist even in central Au-Au collisions 
and dominate (pt) fluctuations in Hijing at lower energies, 
consistent with the jets-off results in the right panel. In 
contrast to Hijing we observe that string-related correla- 
tions in RHIC data are rapidly eliminated with increasing 
Au-Au centrality, even for fairly peripheral collisions 37 1 . 
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VIII. DISCUSSION 

Copious minijet production (a minijet 'plasma' J3J|) has 
been predicted for Au-Au collisions at RHIC [ILI2,llilIl@- 
The theoretical concept of minijets is low-Q 2 partons 
(mainly gluons) produced in the initial stages of relativis- 
tic nuclear collisions Q. The theoretical parton/minijct 
Q range extends from a 'saturation' limit Q s ~ 1 GeV 
for RHIC collisions up to several GeV (semi- hard) 
Minijets are said to carry most of the transverse en- 
ergy in central A-A collisions at RHIC [3] and were 
believed (prior to these measurements) to equilibrate 
rapidly, driving experimentally-observable hydrodynamic 
phenomena (e.g., elliptic flow) Given the theoretical 
uncertainty and expected dominance of minijets in the 
early stages of RHIC collisions it is important to test the 
observability and degree of equilibration of low-Q 2 (Q ~ 
1-5 GeV) partons. 

Is a minijet interpretation allowed for these pt corre- 
lations? A perturbative model of minijet production in 
heavy ion collisions is stated to apply only above par- 
ton p t ~ 2 GeV/c [19(. However, the authors also state 
that the 2 GeV/c lower limit is only a limit on the the- 
oretical description, not the physical phenomenon. And, 
fragments from 2 GeV/c partons should appear at and 
below 1 GeV/c hadron momentum, p-p fragmentation 
functions for parton energies up to 600 GeV measured 
at FNAL extend down to 0.35 GeV/c ^S], and e+- 
e - fra gme ntation functions extend to much lower mo- 
menta [3^!, 40]. The most probable fragment momentum 
in either case is 1-2 GeV/c for a wide range of parton 
energies. Thus, there is no theoretical or observational 
reason which could preclude significant jet fragment con- 
tributions below 2 GeV/c in heavy ion collisions. Con- 
ventional methods for measuring jet angular correlations 
based on a high-p t 'leading particle' are insensitive to 
partons below about 6 GeV. However, the novel analysis 
techniques developed for two-particle correlation analy- 
sis described in this paper and elsewhere [H, [29|, [3(| H3] 
have moved the threshold for direct observation of par- 
tons via final-state hadron correlations down to Q ~ 1 
GeV, and jet-like structure is indeed observed. 

Is a minijet interpretation in fact necessary for these pt 
correlations? The correlation structure in Fig. [2] is domi- 
nated by a same-side peak, an away-side ryA-mdependent 
ridge and a sinusoid. The sinusoid can be interpreted 
as elliptic flow (marking its first observation as a velocity 
phenomenon). A same-side peak (jet cone) and away-side 
ridge are the expected signature angular or number corre- 
lations for high-p t hadron fragments from hard-scattered 
partons. In [7[ it was argued that high-p t angular corre- 
lations on azimuth <p obtained with the leading-particle 
method in A-A collisions are similar to those obtained 
in elementary p-p and e + - e~ collisions with full jet re- 
construction and attributed to hard parton scattering. 
In this analysis the same structure is observed in p t cor- 
relations for p t < 2 GeV/c without a high-p f leading 
particle. Jet structure dominates peripheral collisions in 



Fig. [5] (left panels), and jet-like structure in central col- 
lisions (right panels), although strongly modified, is part 
of a continuous shape evolution from N-N collisions. 

Minijet structures have been observed as angular num- 
ber correlations of low-p t particles in Au-Au collisions at 
130 GeV [13 and p-p collisions at 200 GeV (30(]. Low-p t 
jet-like structures in p t and number correlations are ob- 
served in Hijing Monte Carlo data where the correlation 
mechanism is known to be parton fragmentation [24| . In 
fact, the low-p t jet-like structure that we observe is ex- 
actly what is described as minijets by theory [4l[. We 
conclude therefore that the analogous pt correlations in 
this analysis strongly support a minijet interpretation 
in which hadron fragments from minimum-bias partons 
(no condition is imposed by the analysis on the underly- 
ing parton momentum distribution) are peaked at low 
p t . Ironically, low-Q 2 partons may be more precisely 
and unambiguously characterized by high-statistics pt- 
autocorrelation studies than partons studied in conven- 
tional high-p t leading-particle studies with their biased 
parton momentum spectra and background subtraction 
issues. 

As noted, Ap/^/p re j is within a constant (e r) e^/a^ ~ 
0.25) Pearson's correlation coefficient, a measure of the 
relative covariance between fluctuations in pairs of bins. 
Ap/ ' -s/Pref f° r number correlations measures the num- 
ber of correlated particle pairs per detected particle. If 
(pt) fluctuations are dominated by minijet correlations, 
as our measurements strongly suggest, then increase of 
Ap/ ^/Pref means that either there are relatively more 
minijets (partons) per detected particle with the same 
Pt structure and multiplicity, or the average multiplic- 
ity and/or total pt of minijets has increased, or both. 
That the minijet correlation structure in Fig. [2] is com- 
parable in amplitude to that identified with elliptic flow 
suggests that the two dynamical processes have compa- 
rable importance in heavy ion collisions at RHIC ener- 
gies. The trend with ^/sjvjv in Fig. 2] (left panel) sug- 
gests strong increase of minijet-associated pt production 
at higher energies, and a limit on detectable parton pro- 
duction (correlated hadron fragments) at lower energies, 
with the apparent onset of observable parton fragments 
near y / ijvw = 10 GeV. 

We now return to the question of parton scattering, 
minijets and equilibration in heavy ion collisions. The 
correlation structure of an evolving physical system can 
be used to track the equilibration process. Equilibrated 
or 'thermalized' systems exhibit a large range of corre- 
lation types and degrees, from a Bose condensate to an 
ideal gas. What is relevant for study of a particular sys- 
tem evolving from a non-equilibrium initial state toward 
equilibrium is changes in its correlation structure. In the 
present case parton fragment correlations observed in el- 
ementary N-N collisions should be modified in heavy ion 
collisions depending on the extent and nature of the equi- 
libration process. We can compare observed correlations 
in heavy ion collisions with the expectation for linear su- 
perposition of N-N collisions in a Glauber representation 
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of transparent nuclei as a limiting case. In this analy- 
sis we show that abundant minijet structure associated 
with elementary N-N collisions and exhibiting energy de- 
pendence consistent with QCD expectations for parton 
scattering survives to kinetic decoupling in central Au- 
Au collisions. The structure is strongly modified, but the 
amplitude is still very significant compared to the N-N 
reference. That result suggests that any claims of com- 
plete thermalization for RHIC heavy ion collisions should 
be reconsidered. 



IX. SUMMARY 

In conclusion, we report measurements at several col- 
lision energies of pt angular autocorrelations on pseudo- 
rapidity and azimuth difference variables inferred by in- 
verting the scale dependence of (pt) fluctuations. We 
also report the energy dependence of the centrality and 
scale dependence of (pt) fluctuations to provide compar- 
isons with CERES measurements at SPS energies. The 
Pt autocorrelation distributions are interpreted as con- 
sisting mainly of minimum-bias parton fragments {i.e., 
dominated by minijets). We examine energy dependence 
over the broadest possible interval by comparing STAR 
measurements to compatible CERES measurements at 
the CERN SPS. The results are consistent across ener- 
gies and experiments: 1) excess (pt) fluctuations increase 
fourfold from 20 to 200 GeV for central Au-Au collisions, 
but saturate above 62 GeV for peripheral Au-Au and p-p 
collisions; 2) those fluctuations correspond mainly to pt 



angular correlations identified as minijets, which in cen- 
tral Au-Au collisions are strongly deformed relative to 
p-p collisions (e.g., the same-side peak is broadened on 
^a); 3) (pt) fluctuations thereby associated with initial- 
state parton scattering increase with energy proportional 
to hi{y/sNN} above an onset of detectable parton frag- 
ments near V s aw = 10 GeV. p t autocorrelations thus 
reveal substantial parton fragment correlations surviv- 
ing from initial-state scattering even in central Au-Au 
collisions. Claims of complete thermalization in central 
HI collisions should be reexamined in light of these re- 
sults. The strong energy dependence observed in these 
correlation and fluctuation data should motivate addi- 
tional measurements at lower RHIC energies, spaced in 
energy according to the observed logarithmic trend, to 
investigate the onset of parton fragment production and 
the relationship of low-Q 2 parton scattering to the QCD 
phase boundary. 
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